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ABSTRACT 

We present an analysis of the galaxy distribution surrounding 15 of the most luminous (> 10^* Lq; 
Mi45o ~ —30) QSOs in the sky with z ~ 2.7. Our data are drawn from the Keck Baryonic Structure 
Survey (KBSS), which has been optimized to examine the small-scale interplay between galaxies 
and the intergalactic medium (IGM) during the peak of the galaxy formation era at z ~ 2 — 3. 
In this work, we use the positions and spectroscopic redshifts of 1558 galaxies that lie within '^3' 
(4.2/i~^ comoving Mpc; cMpc) of the hyperluminous QSO (HLQSO) sightline in one of 15 independent 
survey fields, together with new measurements of the HLQSO systemic redshifts. By combining 
the spatial and redshift distributions, we measure the galaxy-HLQSO cross-correlation function, the 
galaxy-galaxy autocorrelation function, and the characteristic scale of galaxy overdensities surrounding 
the sites of exceedingly rare, extremely rapid, black hole accretion. On average, the HLQSOs lie 
within significant galaxy overdensities, characterized by a velocity dispersion cr„ ~ 200 km s~^ and 
a transverse angular scale of ~25" (~200 physical kpc). We argue that such scales are expected for 
small groups with log(Mh/MQ) ~ 13. The galaxy-HLQSO cross-correlation function has a best-fit 

correlation length r^'^ = (7.3 ± 1.3)/i~^ cMpc, while the galaxy autocorrelation measured from the 
spectroscopic galaxy sample in the same fields has r^'^ = (6.0±0.5)/i~^ cMpc. Based on a comparison 
with simulations evaluated at z ~ 2.6, these values imply that a typical galaxy lives in a host halo 
with log(Mh/Mo) = 11.9±0.1, while HLQSOs inhabit host halos of log(Mh/M0) = 12.3±0.5. In spite 
of the extremely large black hole masses implied by their observed luminosities [log(MBH/MQ) > 9.7], 
it appears that HLQSOs do not require environments very different from their much less luminous 
QSO counterparts. Evidently, the exceedingly low space density of HLQSOs (< 10~^ cMpc~^) results 
from a one- in-a- million event on scales << 1 Mpc, and not from being hosted by rare dark matter 
halos. 

Subject headings: galaxies: high redshift - quasars: general - large-scale structure of universe 



1. INTRODUCTION 

The study of galaxies with supermassive black holes 
has become a topic of considerable interest, particularly 
since the discovery that properties of these black holes 
are stron gly correlated wit h those of their host galax- 
ies (e.g. [M agorri an et"n] 119981: IGebhardt erall 120001: 
iFerrarese fc Mcrritt] |200^ The processes of supermas- 
sive black hole accretion and growth can produce spec- 
tacularly luminous QSOs, allowing their study over vast 
cosmological volumes (0 < z < 7). The details of these 
accretion processes, however, are concealed not only by 
distance, but also by our lack of knowledge concerning 
the duty cycle of AGN and the environments that drive 
and sustain their growth. 

Because the brightest QSOs are extreme, ultra- 
luminous objects, it is often assumed that they must in- 
habit comparably rare environments. In particular, the 
rarity of these objects could arise because they require 
the highest mass dark matter halos, which are highly 
biased with respect to the overall matter distribution, 
or because of other, finely-tuned environmental factors 
that influence the availability of gas and the propensity 



^ Based on data obtained at the W.M. Keck Observatory, 
which is operated as a scientific partnership among the Cah- 
fornia Institute of Technology, the University of California, and 
NASA, and was made possible by the generous financial support 
of the W.M. Keck Foundation. 



for the black hole to accrete. As such, the masses and 
spatial distribution of the dark matter halos that host 
QSOs are of considerable interest, and detailed statis- 
tics on these quantities have become available through 
large-scale surveys, primarily through studies of QSO 
clustering using the two-point correlation function. Re- 
cent surveys have covered wide regions of the sky and 
large range s of redshift e.g., the Slo an Digital Sky Sur- 
vey (SDSS; l^rkeFaLlUOOO; Eiscnst eiTetSllmil) . the 
2dF QSO Redshift Survey (2QZ; Groom et al.'*2004f), and 
the DEEP2 Redshift Survey (Davis ct al. 2003). 

Because these surveys include QSOs with a wide range 
of luminosities and redshifts, the QSO autocorrelation 
function has been frequently used to constrain QSO clus- 
tering out_to_Jiigh reds hifts using the SD SS samples 
(e.g. IMvers et al. 20061, [2007t IShen et all [2 007. 201(1 
Ross e t al."2009'). 2QZ samp les (e.g. |P orciani ct al. 2001 



Porciani fc Norbcr g 2006: Groom et al. 2005). and com- 



bined 2dF-SDSS LRG and QSO sur vey (2SLAQ; sur- 
vey description in IGroom et al.l 120091 clustering results 

in Ida Angela et ai1l2008D . The results of these analyses 
are in broad agreement that QSOs inhabit host dark mat- 
ter halos of mass log(Mh/M0)~12.5 at redshifts z < ?,. 
Due to the low space density of QSOs at all redshifts, 
these autocorrelation measurements have generally been 
confined to large scales, but co mplementary m e asure - 
ments have also been obtained: IHennawi et al.l ()2006f ) 
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and iShen et all ()201Cl[ ) conducted surveys for close QSO 
pairs; the galaxy-QSO cross-correlatio n function was 
measured bv lAdelberger fc S teidel' C2005[). in the DEEP2 
survey at z ~ 1 bv lCoil et al. (2007), a nd in a low redshift 
z < 0.6) SDSS QSO sample by Padmanabha n et alJ 
20091 ). These studies generally agree with the QSO auto- 
correlation results, and the mass scale log(Mh/MQ)'^12.5 
seems fairly well-established for the general population of 
QSOs at z < 3. 

However, studies which divide the population of QSOs 
into specific subsamples reveal a more complicated pic- 
ture of the dependence of QSO properties on halo mass. 
Low-redshift studies display a possi ble relation betwee n 
obscuration and host halo mass (|Hickox et all [20TTI) . 
which may be significant at higher redshifts, where the 
population of obsc ured QSOs is relatively unconstrained. 
IShen et al.l (poToh find that radio-loud QSOs are more 
strongly clustered than radio-quiet QSOs matched in red- 
shift and optical luminosity. In addition, there is an 
expected dependence of QSO luminosity on host halo 
mass because the QSO luminosity depends on black hole 
mass, which in turn exhibits the aforementioned associ- 
ation with the mass of the host halo. In practice, how- 
ever, QSO luminosities depend in detail on the availabil- 
ity of matter to accrete and the physical processes gov- 
erning the efficiency with which this accretion occurs. 
Thus, it is perhaps not surprising that the clustering of 
QSOs shows little associati on with QSO luminosity in 
observations near z ~ 2 (e.g. lAdelberger fc Steidell[2nn5l : 

ICroom et al.l 120051: Ida Angela et al.ll2008l and in siru ula- 
tionsbv lLidz et al.ll2006[ ): however. IShen et al.l (|2010[ ) de- 
tect stronger clustering among the most lumi nous QSOs 
in their sample at z > 2.9, and Kru mpe et all (2010) find 
that SDSS QSOs at z ^ 0.25 cluster more strongly with 
increasing X- ray luminosity. Finally , the survey of close 
QSO pairs bv lHennawi et al.l (|2006l ) reveals an excess at 
the smallest scales, which the authors attribute to dis- 
sipative interaction events that trigger QSO activity in 
rich environments. In short, the properties of QSOs are 
related to their host halo masses in a complex manner, 
and it is clear that other environmental factors are in 
play. 

In this paper we study the environments of hyperlu- 
minous QSOs (HLQSOs; defined here by a luminosity 
log(zyL,y/L0)>14 at a rest-frame wavelength of 1450 A) 
at 2.5 < 2 < 3 by measuring the magnitude and scale 
of overdensities in the galaxy distribution at small (^3') 
projected distances using data from the Keck Baryonic 
Structure Survey (KBSS). This approach complements 
existing studies in numerous ways: targeting narrow 
fields allows us to study the local environments of these 
extremely rare HLQSOs, including the galaxies at com- 
parable redshifts that lie far below the flux limits of the 
typical wide-field QSO surveys. In this way we are able 
to constrain the properties of the relatively unexplored 
environments of the highest-luminosity QSOs. Focus- 
ing on the brightest QSOs should reveal whether host 
halo mass plays a significant role in determining QSO 
properties, while sensitivity to the local environment may 
demonstrate whether these HLQSOs are associated with 
the types of environments where mergers and dissipative 
interaction are expected to be most common. 

This paper is organized as follows: in Sj2]we discuss the 



observations used in this study; f^l describes the tech- 
niques used to construct an unbiased measure of the 
galaxy distribution around the HLQSOs and our esti- 
mates of the magnitude and scale of the surrounding 
galaxy overdensities. In ^ we describe and implement 
a method for estimating the small-scale galaxy-HLQSO 
correlation function and galaxy-galaxy autocorrelation 
function from our data along with the implied galaxy 
and HLQSO host halo masses. In ij5]we present evidence 
that the HLQSOs inhabit group-sized virialized struc- 
tures conducive to merger events; a summary is given in 
f|6l Throughout this paper, we will assume rim = 0.3, 
= 0.7, and h = i/o/(100 km s'^). We have left ah 
comoving length scales in terms of h for ease of compar- 
ison to previous studies, but we quote physical scales, 
luminosities, and halo masses assuming h = 0.7. For fur- 
ther clarity, we denote comoving distance scales in units 
of cMpc (comoving Mpc) and physical scales as pkpc 
(physical kpc). 

2. DATA 

The data used in this study form part of the Keck 
Baryonic Structure Survey (KBSS; Steidel et al. 2012), a 
large sample (Ngai — 2298) of high-redshift star-forming 
galaxies (1.5 < z < 3.6) close to the hnes-of-sight of 15 
HLQSOs at redshifts 2.5<z<2.9. Because we have ob- 
served fields of differing solid angle around each of these 
HLQSOs, we standardized the fields for the purposes of 
this study by including only those galaxies within 66 ^ 3' 
(4.2/1^1 cMpc at the HLQSO redshifts) of the hnc-of- 
sight of the HLQSO in each, an area that is well-sampled 
for all 15 fields. This subset of the total KBSS dataset 
contains 1558 galaxies and comprises the entire sample 
used in this paper. 

2.1. HLQSO Redshifts 

An important prerequisite to establishing the galaxy 
environment of the HLQSOs is an accurate measure- 
ment of the HLQSO systemic redshifts. Redshifts for 
QSOs in the range 2 < ^^qso ^ 3 are typically mea- 
sured from the peaks or centroids of broad emission lines 
of relatively high ionization species in the rest-frame 
far-UV (e.g., N V A1240, C IV A1549, Si IV A1399, 
C III] A1909). These lines are known to yield redshifts 
that differ significantly from systemic, and tend to be 
blue-shifted by several hundred to several thousand km 
s~^ (see e.g. iMcIntosh e t al. 1999; Richards ct al. 2002!; 
lGon(;alves et al.l2008f ). These velocity offsets also tend to 
increase with QSO luminosity, thus making the present 
sample of hyper luminous QSOs particularly susceptible 
to this issue. In view of the importance of precise red- 
shifts to locate the HLQSO environments within the sur- 
vey volume, we obtained near-IR spectra of the entire 
sample using NIRSPEC on the Keck II 10m telescope, 
TripleSpec on the Palomar 200-inch (5m) telescope, and 
in some cases, both (see Table [T]). 

Among the 15 HLQSOs in the sample, narrow forbid- 
den lines ( [O III] A5007) were detected for only 2 of them 
(Q1623-t-268, Q2343+12), either because no such lines 
were present in the spectra (common at the highest lumi- 
nosities), or because the HLQSO redshift was such that 
the strongest transitions fell in regions between the near- 
IR atmospheric bands. However, in all cases we were able 
to measure one or more hydrogen Balmer lines and the 
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TABLE 1 

HLQSO Redshifts and Corrections 



QSO NIR Spectra Source'' ^ncw'' ^old" Av (km 8"^) 



QOlOO+13 (PHL957) 


Keck II/NIRSPEC 


2.721 


±0, 


,003 


2, 


,681 


-0, 


,040 


-3214 


HS0105+1619 


P200/TSPEC 


2.652 


±0, 


,003 


2, 


,640 


-0, 


,012 


-983 


Q0142-10 (UM673a) 


Keck II/NIRSPEC 


2.743 


±0, 


,003 


2, 


,731 


-0, 


,012 


-943 


Q0207-003 (UM402) 


P200/TSPEC 


2.872 


±0, 


,003 


2, 


,850 


-0, 


,022 


-1699 


Q0449-1645 


P200/TSPEC 


2.684 


±0, 


,003 


2, 


,600 


-0, 


,084 


-6818 


Q0821+3107 (NVSS) 


P200/TSPEC 


2.616 


±0, 


,003 


2, 


,624 


+0, 


,008 


+686 


Q1009+29 (CSO 38) 


Keck II/NIRSPEC 


2.652 


±0, 


,003 


2, 


,620 


-0, 


,032 


-2620 


SBS1217+490 


P200/TSPEC 


2.704 


±0, 


,003 


2, 


,698 


-0, 


,006 


-484 


HS1442+2931 


P200/TSPEC 


2.660 


±0, 


,003 


2, 


.638 


-0, 


,022 


-1797 


HS1549+1919 


Keck II/NIRSPEC 


2.843 


±0, 


,003 


2, 


.830 


-0, 


,013 


-1011 


HS1603+3820 


P200/TSPEC 


2.551 


±0, 


,003 


2, 


.510 


-0, 


,041 


-3452 


Q1623+268 (KP77)'^ 


Keck II/NIRSPEC 


2.5353 


±0, 


,0005 


2, 


.518 


-0, 


,018 


-1489 


HS1700+6416 


Keck II/NIRSPEC 


2.751 


±0, 


,003 


2, 


,736 


-0, 


,015 


-1220 


Q2206-199 (LBQS) 


Keck II/NIRSPEC 


2.573 


±0, 


,003 


2, 


,558 


-0, 


,015 


-1255 


Q2343+12 (also SDSS) 


Keck II/NIRSPEC 


2.573 


±0, 


,003 


2, 


,515 


-0, 


,058 


-4854 



Refers to the instrument used to measure the near-IR QSO spectra and redshift. NIRSPEC is 
used on the Keck II telescope, while P200 is the Palomar Hale 200-inch telescope, used with the 
TripleSpec instrument. 

^ Zncw refers to the redshift used in this analysis. 
2old refers to the previous published redshift value. 

The redshift for Q1623+268 (KP77) is more tightly constrained because of the presence of narrow 
[O III] lines at the presumed systemic redshift of the QSO. 



TABLE 2 

Galaxy Samples and HLQSO Properties 



Field 


^QSO" 


1^1450 


Mbh^ 
(IO^Mq) 






^CDM 


Ntot 


^ 1500^^ 


QOlOO+13 


2.721 


6.4 


2.0 


68 


12 


15 


95 


7 


HS0105+1619 


2.652 


4.5 


1.4 


74 


6 


23 


103 


7 


Q0142-10'= 


2.743 


<6.4 


<2.0 


75 


13 


16 


104 


1 


Q0207-003 


2.872 


6.1 


1.9 


54 


12 


27 


93 


7 


Q0449-1645 


2.684 


4.0 


1.3 


68 


12 


31 


111 


9 


Q0821-|-3107f 


2.616 


4.1 


1.3 


64 


7 


21 


92 


4 


Q1009-I-29 


2.652 


10.9 


3.4 


54 


19 


43 


116 


8 


SBS1217-(-490 


2.704 


5.1 


1.6 


67 


14 


11 


92 


3 


Q1442 -1-2931 


2.660 


4.9 


1.5 


71 


25 


22 


118 


3 


HS1549+1919 


2.843 


14.9 


4.6 


54 


14 


39 


107 


23 


HS1603+3820 


2.551 


11.0 


3.4 


80 


15 


14 


109 


10 


Q1623-KP77* 


2.5353 


3.2 


1.0 


82 


9 


12 


103 


7 


HS1700-(-6416 


2.751 


13.6 


4.3 


69 


16 


16 


101 


6 


Q2206-199 


2.573 


4.5 


1.4 


78 


11 


20 


109 





Q2343-I-12 


2.573 


3.8 


1.2 


71 


9 


25 


105 


6 



^ 2Qgo refers to the systemic redshift of the field defined by the HLQSO (see Table [T] and 

^ Li45o refers to the estimated luminosity uLi, near a rest-frame wavelength X^cst — 1450A, 
extrapolated from the g' and r' magnitudes from the SDSS Q^isenstein et al. 2011) database 
when available, and otherwise from our own measurements. We have assumed h = 0.7. 

Mgjj is the minimum black hole mas s cap able of producing a QSO with luminosity L]^45o, 
assuming Eddington-limited accretion ( t|4.6ll . 

^1500 ~ N(|<5t;| < 1500 km s~^) is the number of galaxies in the field that have spectroscopic 
redshifts within 1500 km s~^ of their corresponding HLQSO. 

Q0142-10 (UM673a) is known to be gravitationally lenscd (Surd ei et al.l [19871 ) and has an 
unknown magnification; the estimated luminosity and mass arc therefore upper limits. 
* Q0821-I-3107 and Q1623-KP77 arc the onl y HLQSOs in our sam ple with radio detections. 
Q0821 has a fiux = 162 mjy at 48 30 MHz IILangston et al.|[T990l) : KP77 has a flux /„ = 6.4 
mjy at 1.4 GHz l lCondon et aI.|[T99a) . 
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Mg II A2798 line, which were mutually consistent and are 
known to be closer to the true syst emic redshift than the 
high ionization lines in the UV ()McIntosh et al.l 119991 : 
[Richards et alJ [20021 ) . The redshifts obtained from the 
Balmer/Mg II lines (which agree well with that given 
by [O III] in the two cases where all were measured) 
were then subjected to several cross-checks, including: 
the wavelength at the onset of the Lyman-a forest mea- 
sured from the high resolution HLQSO spectra (a lower 
limit on the systemic redshift, but one which agrees to 
within Az ~ 0.001 of the Balmer line redshift in all but 
2 cases); the redshift of narrow He II A1640 in inter- 
mediate resolution optical spectra of the HLQSOs; and, 
in several cases, regions exhibiting narrow Lya emission 
were discovered with small angular separations from the 
HLQSO, and we have found that such nebulae lie very 
close to the systemic redshift of the nearby HLQSO. In 
2 cases (HS1603+3820 and Q1009+29) this last crite- 
rion led to a significant modification (Az ^ -fO.Ol, or 
~ 800 km s^^) of the redshift suggested by the near-IR 
spectroscopy. We adopt a HLQSO redshift uncertainty 
(Tz = 270 km s~^ (for those HLQSOs without measured 
[O III] redshifts) based on the measure d dispersion of 
the M g II line with respect to [O III] by [Richards et al.l 
(|200l : the broad agreement among our many redshift 
criteria suggest that this is a conservative estimate of 
the redshift uncertainties. 

Table [T] summarizes the adopted redshifts for all 15 
HLQSOs based on these considerations; also given (col- 
umn 4; Zoid) is the published redshift for each and the 
redshift and velocity error that would result from adopt- 
ing the published values (Az = Zou—Znew)- As expected, 
all but one of the old redshifts are systematically too low 
(the median shift is ~ —1500 km s^^, and the mean 
~ —2100 km s^"'^). Failure to account for these large 
velocity errors would severely compromise our measure- 
ments. As we show below, the measured zqso values 
must be quite accurate given the very tight redshift-space 
correlation between the HLQSOs and the spectroscopi- 
cally measured, continuum-selected galaxies nearby. 

2.2. Galaxy Redshifts 

Galaxy redshifts were measured using low-resolution 
(~5A), rest-frame UV spectra obtained with the LRIS 
multi-object spectrograph on the Keck I telescope 
(|Oke et al.l [19951: iSteidel et al.l 120041) . Candidate galax- 
ies were color-selected using the Lyman-break technique 
and were sorted as BX (z 2.2), MD (z ~ 2.6), or CDM 
(z ^ 3) galaxies base d o n the color criteria disc ussed 
in ISteidel et"all (|2003f) and lAdelberger et al.l ()2004D : the 
data collection and reduction procedures are described 
therein. All galaxies in the spectroscopic sample have 
TZ < 25.5 [where TZ = mAB(6830A)], which corresponds 
to Mab(1700A) < -19.9 at z - 2.7 [about 1 magnitud e 
fainter than at this redshift; see iReddv et al.l [20081) . 
Redshifts were determined by a combination of Lya emis- 
sion or absorption and far-UV interstellar (IS) absorp- 
tion. Since Lya emission tends to be redshifted with 
respect to the systemic redshift of the host galaxy, and 
interstellar absorpt i on tends to be bluesh i fted (see e.g. 
Shaplev et alll2003l: lAdelberger et al.ll2003HSteidel et al] 



20101) . we estimate each galax y's systemic redshift via 
the method proposed in lAdelberger et al.l ((2005aj) and 



updated by ISteidel et al.l ()2010D . In this method, the 
average Lya emission and IS absorption offsets are cal- 
culated based on the redshift of the Ha nebular line (NIR 
spectroscopy is available for a subset of the galaxy sam- 
ple), which traces ionized gas in star- forming regions of 
the galaxy, and is thus a more accur ate estimate of the 
systemic redshift. iRakic et al.l (|2011[ ) derive similar cor- 
rections for the same galaxy sample using the expected 
symmetry of IGM absorption about the systemic redshift 
of the galaxy. 

We estimate the systemic galaxy redshifts (zgai) based 
on a combination of the above results. A more detailed 
discussion of our correction formulae can be found in 
Rudie et al. (2011; in prep.), but the formulae are re- 
produced below. For galaxies with NIR spectra (e.g. the 
Ha line), the NIR redshift is used with no correction. 
For galaxies with measured Lya emission but without 
interstellar absorption, we use the following estimate: 



(1) 



where ZLya is the redshift of the measured Lya emission 
and Avhya — —300 km s~^ is the velocity shift needed 
to transform the Lya redshift to the systemic value, Zgai. 

For galaxies with interstellar absorption, we use an es- 
timate based on the absorption redshift whether or not 
Lya emission is present: 



Zgai = Zis 



Avis 



(l + ^is) 



(2) 



where zis is the redshift of the measured interstellar ab- 
sorption and AwLya = 160 km s^^ is the velocity shift 
needed to transform the absorption redshift to the sys- 
temic value. 

For galaxies with both interstellar absorption and Lya 
emission, we verify that the corrected absorption redshift 
does not exceed the measured redshift of the Lya line; 
that is, we verify that zjg < Zgai < ZLya, where Zgai is 
calculated using Eq. [5] above. If this condition is not 
satisfied, we recompute the galaxy systemic redshift as 
the average of the absorption and emission redshifts: 



ZlS + ^Lya 



(3) 



The residual redshift errors (calculated from the galax- 
ies in the NIR sample) have a standard deviation 
fv.crr — 125 km s^^, which we adopt as the uncertainty 
in our galaxy redshift measurements. 

3. REDSHIFT OVERDENSITY 

In order to consider the positions of the galaxies rel- 
ative to their corresponding HLQSOs in redshift space 
while accounting for the differences in the HLQSO red- 
shifts between fields, the redshift of each galaxy was 
transformed into a velocity relative to its associated 
HLQSO. For a galaxy with index i in a field with index 
J, this velocity difference is given by 



5vi 



1 



(Zgal,! — ZqsOj) • 



(4) 



Once transformed to units of velocity, the distributions 
of galaxies relative to their HLQSOs were stacked to re- 
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Fig. 1. — Velocity distribution of galaxies with respect to their 
nearest HLQSOs, stacked for all 15 HLQSO fields. The velocity 5v 
is given by Eq. |4] where 5v = for a galaxy at the redshift of its 
corresponding HLQSO. The yellow shaded area c orre sponds to the 
selection function, constructed as described in i|3.1l The dashed 
curve is a gaussian profile fit to the overdensity, with cr^ gj. = 350 

km s~^. After removing the effect of our ctv err ~ 125 kms"! (270 
km s~^) galaxy (IfLQSO) redshift errors, we estimate a peculiar 
velocity scale of (Tv,pec — 200 km for the galaxies associated 
with the overdensity, with an offset (Sv) = 106±54 km s"'^ from 
the HLQSO redshifts. 

veal the average environment of HLQSOs in terms of the 
local galaxy number density (per unit velocity)~this dis- 
tribution is shown in Fig. [TJ The distribution shows a 
well-defined peak near Sv = 0, indicating the presence of 
significant clustering of the galaxies around the HLQSO 
redshifts. We attribute the slight offset of the overden- 
sity from the HLQSO redshifts (fit (6v) = 106±54 km 
s~^) to a residual systematic offset in our determination 
of the HLQSO redshifts. 

3.1. Building the Selection Function 

Clustering measurements can be grossly misinterpreted 
when the r elevant selection functions are not well- 
understood (|Adelbergeiil2005l) . While the criteria for se- 
lecting galaxies for follow-up spectroscopy were identical 
for all 15 of the KBSS fields, small differences in image 
depth and seeing, as well as slight changes in the algo- 
rithms used for assigning relative weights in the process 
of designing slit masks, can lead to field-to-field vari- 
ations in the redshift selection functions. To at least 
partially mitigate such variations in the redshift-space 
sampling between fields, we used the number of success- 
fully observed BX, MD and CDM galaxies in each field 
to estimate the form of our field-specific selection func- 
tions J^j{z). These estimates of the selection functions 
were constructed as follows. 

First, the redshift distributions of all BX, MD and 
CDM galaxies in our sample were arranged in a coarse 
histogram with bins of width Az = 0.2. A spline fit 




Redshift 

Fig. 2. — Redshift distributions for BX, MD, and CDM color- 
selected galaxy types. Red histograms display the measured distri- 
butions of all such galaxies in each sample, while the yellow region 
represents the fit spline function specific to the color- selected sam- 
ple (ie. Af-Qxi^)' ■^MDi^)^ ^''^'^■^CDMi^))- The blue hashed region 
is the overall redshift selection function for all color types. 

was then performed to estimate the smooth distribution 
functions of each galaxy type-the histograms and spline 
fits for each type are displayed in Fig. [21 

For each field, we built a field-specific selection function 
by combining these galaxy redshift distributions for each 
color criterion according to the number of those galaxies 
successfully observed in the field. Thus for a field with 
index j, the redshift selection function is given by Eq. [5j 

NMI),j■^fMu{z) + 

NcDM,3^fcDM{z) (5) 

where A^BX,j corresponds to the number of BX-selected 
galaxies in field j, TVbx is the selection function for BX- 
selected galaxies over all fields, and other variables are 
defined similarly. We then transform these redshift-space 
selection functions into units of velocity relative to their 
corresponding bright HLQSOs using Eq. |H Finally, we 
combined this set of field-specific velocity-space selection 
functions (already weighted by the number of galaxies in 
each field) into a single stacked function: 

mv)^T.^jiv)- (6) 

The resulting selection function is fairly flat over the 
range \Sv\ < 20000 km s^^ with a slight negative slope 
(yellow shading in Fig. [IJ, indicating our slight bias to- 
ward detecting objects "in front" of the HLQSO (that 
is, at lower redshifts) compared to galaxies slightly "be- 
hind" the HLQSO in each field. The selection function is 
thus a prediction for the observed distribution of galaxies 
in relative- velocity space in the absence of clustering. 
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3.2. Bias in Field Selection 

We previously knew one KBSS field (HS1549+1919) 
to have a large overdensity in the galaxy distribution 
very close to the redshift of the central HLQSO. The 
variation in overdensity among fields can be estimated by 
-^1500 in Table [2l which is the number of galaxies within 
1500 km s-i of the HLQSO redshift for that field. In 
order to ensure that our clustering results are not being 
dominated by a single field, we repeated our analysis on 
subsamples of the data consisting of 14 of the 15 fields, 
removing a different field each time. In each case the 
magnitude and scale of the overdensity was consistent 
with that observed when all 15 fields were included in 
the analysis, indicating that the observed magnitude and 
scale of the overdensity arc not determined by any single 
field. 

Hubble Distance from QSO (h~' comoving Mpc) 
-100 100 



/ovr — (-^obs — ^prcd)/^prcd 



(7) 




6v (10" km s"') 

Fig. 3. — The relative overdensity /ovr (Eq. O as a function of 
velocity relative to the central HLQSOs over a wide velocity range. 
The overdensity is measured in bins of 500 km s"'^, as this Hubble 
flow velocity roughly corresponds to the same physical scal e as o ur 
transverse field of view {5h~^ cMpc ~ 500 km s~^). See ^3.11 for 
details on the selection function. 



3.3. Redshift Clustering Results 

Fig. [T] shows the observed galaxy distribution in units 
of velocity along with the selection function estimate 
from i j3.1l The peak in the galaxy distribution near the 
HLQSO redshifts is clearly visible. Fitting a Gaussian 
function to the histogram in Fig. [l]gives a velocity width 
Cv,fit = 350±50 km s~^, which includes the effect our 
(Tv,err — 125 km s"^ galaxy redshift errors and the ran- 
dom residual errors in our HLQSO redshifts, assumed to 
be Cv.err ~ 270 km s ^. After subtracting the redshift 
errors in quadrature, we find an intrinsic velocity width 
of (Tv.pcc — 200 km s^^ for the galaxy overdensity, which 
we attribute to peculiar velocities. Note that the resid- 
ual HLQSO redshift errors are uncertain and likely to be 
largely systematic (see ^2.1\i . so our estimated velocity 
dispersion is an upper limit on the true peculiar velocity 
scale if the random component of the HLQSO redshift 
error is larger than we have assumed. 

We also consider the relative overdensity at the 
HLQSO redshift by comparing the observed density to 
that predicted by our selection function. The distribu- 
tion is plotted as a relative overdensity 



in Fig. [21 where iVobs is the number of galaxies observed 
in a given velocity bin and A^piod is the number predicted 
for that bin by our selection function. The relative over- 
density is measured in bins with Av = 500 km s~^; this 
scale was chosen to correspond roughly to the transverse 
scale of our field, since a Hubble-fiow velocity of 500 km 
s^^ ^ 5h~^ cMpc at these redshifts. Fig. [3] shows that 
the HLQSOs are associated (on average) with a Sn/n ^ 7 
overdensity of galaxies when considered on the ~5ft.~^ 
Mpc scale of our field, with no features of comparable 
amplitude over a wide range of redshifts (40000 km s~^ 
corresponds to Az ~ 0.5 at z ~ 2.7). 

Angular Distance from QSO (') 
1.0 1.5 2.0 2.5 




12 3 4 

Projected Distance from QSO (h"' cMpc) 

Fig. 4. — For each projected circular annulus, the relative over- 
density {/ovr; Eq. [7]l of galaxies within 1500 km of the HLQSO 
with respect to the redshift selection function and angular selection 
function. The overdensity of galaxies is localized for the most part 
to a tranverse scale R < 0.5h^^ cMpc. 

Repeating this analysis after dividing the galaxies into 
radial annuli, we find that the redshift association is most 
pronounced for those galaxies within 25" of the HLQSO 
line of sight ('--^200 pkpc), though a lower level of red- 
shift clustering does extend to larger projected distances 
(see Fig. |4|). If this distance is taken as an isotropic 
spatial scale of the galaxy overdensity, then the line-of- 
sight velocity dispersion due to the Hubble fiow would 
be only ~65 km s"^. However, a less-significant overden- 
sity does extend to larger radii, and thus likely includes 
many galaxies that are clustered around the HLQSO but 
move with the Hubble fiow. In order to ensure that the 
measured velocity width is not infiated by these non- 
virialized galaxies, we directly measure the velocity dis- 
persion among the 15 galaxies within 1500 km s~^ and 
0.5 h^^ cMpc (200 pkpc) of the HLQSOs; as discussed 
in ij5l these galaxies are likely to be virialized and as- 
sociated with the HLQSO, and our selection functions 
predict only 1.5 galaxies in this volume in the absence 
of clustering. These 15 galaxy velocities have a sample 
standard deviation of 335 km s~^, consistent with the 
velocity width measured for the entire overdensity. The 
observed velocity spread is thus presumably set by pecu- 
liar velocities of (Tv.pcc — 200 km s"^ among the HLQSO- 
associated galaxies. 
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A comparison of Figs. [3] & H] demonstrates that the 
relative overdensity is highly scale-dependent. If we as- 
sume that the width of the overdensity in velocity space 
is entirely due to peculiar velocities, and hence that all 
15 of the galaxies observed with R < 0.5h~^ cMpc and 
< 500 km s~^ are physically located within a three- 
dimensional distance r < 0.5/i~^ cMpc from their nearest 
HLQSO, then the number of galaxies in this composite 
volume is ^5Gx the number predicted by our redshift and 
angular selection functions (described in i j3.1l & iJ4.1.11 
respectively) . 



P{R)dn = P^dn[l + Wp{R)]=dn [ P{r)dz (8) 

where R — Da{z)9{1 + z) in the projected comoving dis- 
tance from the HLQSO, and Da{z) is the angular diam- 
eter distance to the HLQSO. In the limit Az oo, and 
assuming a power-law form of the three-dimensional cor- 
relation function, it can be shown that the reduced angu- 
lar correlation function has an equally simple power-law 
form; 



Wp{R) = AR-"^ 



(9) 



4. CORRELATION FUNCTION ESTIMATES 

4.1. Galaxy-HLQSO Cross- Correlation Function 

Much of the recent work on QSO clustering relies on 
large-scale two-point correlation functions, parti cularly 
the Q SO autocorrelation function (see e.g. Shen et ahl 
\mrj\\ . The galaxy-HLQSO cross-correlation function 
can provide a complementary estimate of HLQSO host 
halo mass. 

The correlation function is defined as the excess condi- 
tional probability of finding a galaxy in a volume dV at 
a distance r = |ri — r2|, given that there is a HLQSO at 
point ri, such that P{r2\ri)dV = Po[l + £,{r)]dV, where 
PodV is the probability of finding a galaxy at an average 
place in the universe. Here we assume a power-law form 
for the correlation function: ^gq = j where 7 

is the slope parameter and corresponds to the comov- 
ing distance at which the local number density of galaxies 
is twice that of an average place in the universe. 

Many recent analyses of the two-point correlation 
function have dispensed with power-law fits in fa- 
vor of directly model ing the halo -occupatio n distribu- 
tion ( HOD; see e.g . iSeliakI [2000t [Berlind fc Weinberd 
2002t iZehavi et al.l I2004D based on the theory of 



Press fc Schechtenlj 19741 ) and a statistical method of 
populating dark matter (DM) halos with galaxies. A 
general feature of these HOD models is a deviation from 
a single power law at distances near lh~^ cMpc due to a 
transition from the single-halo regime (the clustering of 
galaxies/QSOs within a single dark-matter halo) to the 
two halo regime (the clustering of galaxies/QSOs hosted 
by distinct halos). 

In this paper we implement the simpler power-law 
fitting technique for the following reasons. First, our 
smaller sample (with respect to the large surveys at low 
redshift) does not allow us to detect a deviation from a 
power-law fit with any significance, particularly for the 
galaxy-HLQSO cross-correlation. Second, our choice to 
fix the power-law slope 7 (see below) desensitizes our 
result to the precise shape of the correlation function, 
leaving the clustering length r^'^ to primarily reflect the 
integrated pair-probability excess over the range of pro- 
jected distances in our sample. 

In practice, the three-dimensional correlation function 
^(r) is not directly measurable: line-of-sight velocities 
are an imperfect proxy for radial distance due to pecu- 
liar velocities and redshift errors. As such, it is more 
useful to consider the reduced angular correlation func- 
tion, Wp{R\Az) by integrating over a redshift or velocity 
window: 



with 77 = 7 — 1. However, we would like to restrict 
our analysis to small redshift /velocity scales, choosing 
a value of Az that includes the entire clustering sig- 
nal while eliminating the noise contribution of uncorre- 
lated structure at large line-of-sight separations from the 
HLQSO, and this priority precludes the assumption of 
Az — 00. In the case of a truncated redshift range, the 
reduced angular correlation function does not simplify to 
a power-law, and instead takes the form of a Gaussian 
hypergeometric function, denoted as 2-^1(0, fo; c; z). In 
our particular case, the reduced angular cross-correlation 
function w^'^ is expressed by the following: 



w^^iR)- 




(VR'^ + z^/r^^ydz 



(10) 



where zq is the half-width of the redshift window over 
which we compute the clustering strength. We choose 
a value zq = 1500 km ~ 14/i~^ cMpc in order 
to encompass the entire observed overdensity (see Fig. 
[1]) and the range of projected distances we are able 
to probe {R < A.2h~^ cMpc) without including excess 
noise. We then fit the reduced angular correlation func- 
tion Wp^{R\r^'^,^) to the data by variation of the cor- 
relation length r^'^ . We fix 7 = 1.5 for simplicity in 
matching our data to halo populations (see §4.3p : this 
value of 7 was chosen as it is a reasonably good fit to both 
the galaxy autocorrelation function and galaxy-HLQSO 
cross-correlation function, as well as the correlation func- 
tions among the simulated halo populations. Increasing 
the value of 7 causes the best-fit value of r^^ to decrease, 
but the corresponding halo mass is very insensitive to the 
choice of 7, so long as the same value is used for both 
the galaxy and the simulated halo populations. 

In order to estimate Wp'^{R) from our data, we sep- 
arate our fields into projected circular annuli of varying 
widths, constructed so that each annulus has a roughly 
similar signal-to-noise ratio, and with our largest an- 
nulus having its outer edge ~200" from the HLQSO, 
a projected distance of i? = A.2h~^ cMpc. As noted 
above, we wish to restrict our analysis to those galax- 
ies that are closely associated with a HLQSO in red- 
shift as well as projected position, so we also separate 
our galaxy sample into two velocity groups: one with 
\6v\ < 1500 km s~^ and one with \dv\ > 1500 km s~^ 
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In this way, we define Ny{Rk) as the number of velocity- 
associated galaxies in the k**^ annular bin and A^o(^fe) as 
the number of non-associated galaxies in the binQ The 
solid angle covered by the k"^ annular bin is given by 
Ak = 7r(i?^^^^j, ]^ - i?f„ncr,k); thus we likewise define the 
area densities of galaxies I]„_o(^fc) = ^v,o{Rk)/Ak- 

We then used the selection function constructed in 
N3.1I to estimate the expected number of galaxies in 
each velocity group [ie. Afy = A/'(|5w| < 1500) and 
A/o = A/'(|(5w| > 1500)], which we convert to expected 
average area densities for each velocity group (S^'^^'^ = 

A/'t),o/-4fioid)- Finally, we divide the measured area den- 
sities by the average predicted value to define the rel- 
ative overdensity of each annulus. If the overdensity 
^v{Rk) /^v'^'"^ is purely due to the clustering signal, then 
the reduced angular cross-correlation function is given 
by u;GQ = S^S,s(-Rfc)/S^^,,ed - 1; however, this assump- 
tion is invalid if the angular sampling of the field is not 
uniform, which we explore below. 

4.1.1. Angular Selection Function 

The use of the redshift selection function in A^prcd en- 
sures that large-scale variations and sampling biases in 
redshift are taken into account in our analysis. Selec- 
tion biases can also occur in the plane of the sky; be- 
cause our fields are centered on their HLQSOs, any bias 
that varies with distance from the center of the field will 
mimic a change in the correlation function. To account 
for this effect, we recall that galaxies with \Sv\ > 1500 
km show no association with the HLQSO (see Fig. 
[T]), and therefore should be uniformly distributed on av- 
erage. Therefore, if the function So(^) is not a con- 
stant, it must describe a non- uniform angular selection 
function, which encapsulates variations in optical selec- 
tion sensitivity (e.g. due to non-uniform extinction or 
field coverage) as well as any biases in slit positions on 
our masks. We assume these biases are independent of 
redshift, and that they produce the same fractional ex- 
cess of galaxy counts in all velocity bins. Therefore, the 
measured values of St,(i?)/SP'^°'^ correspond to the the 
true reduced correlation function 1 -I- Wp^R) multiplied 
by a transverse (angular) selection function, which we 
estimate by T,Q{Rk). 



T,i,{Rk) So(i?) 



-^prcd 



[1 + WpiRk)] 



(11) 



We found that the measured values of Eo(i?fc) are well- 
matched by a power-law in R, and therefore, rather than 
using Eq. [TT] directly to estimate Wp{Rk), we found best- 
fit parameters a and /3 for the following model: 



Eo(i?) 



->pred 



aR^ 



(12) 



The best-fit parameters for this model are a = 1.59, 
(3 = —0.58 with R in h^^ cMpc; the fit selection func- 
tion is displayed in Fig. [S] Combining Eqs. [TUl [TTl & 
1121 we arrive at an explicit model for Y,y(R) in terms of 

^ In this paper, we will use the subscript or superscript v to 
denote velocity-associated galaxies, and to denote non-associated 
galaxies. 



the galaxy-HLQSO cross-correlation length tq and cor- 
relation slope 7: 



aR" 



1 



2F1 



3 

2' 



,2 \ 1 



i?2 



(13) 



where zq = (1500 km s^^)Hq ^l + z)"^ is the half-width 
of the redshift window in physical units, and again a and 
P are set by fitting So(i?) to So(i?fe). We then adjust the 
free parameter tq corresponding to the cross-correlation 
function to fit the measured values of Si,(i?fe). 

The fit to was performed via a simple x^- 
minimization using an error vector constructed assuming 
Poisson uncertainties in the galaxy counts. The binned 
data cover a range of projected distances 0.22 — 3.57h~^ 
cMpc. 

Our empirical estimate for Wp{R) obtained via the 
above methods is displayed in Fig. \6\ We find a best-fit 
correlation length rg = (7.3 ± 1.3)h~^ Mpc after fixing 
7 = 1.5, where the error is a la uncertainty computed 
via a bootstrap estimate. This procedure consisted of 
repeating the entire analysis 100 times (computation of 
selection functions, counting of pairs, and x^-fitting of 
Wp) using a random bootstrap sample of 15 of the 15 in- 
dependent fields selected with replacement. The quoted 
uncertainty is the standard deviation of parameter values 
derived from these 100 bootstrap samples. The results of 
this procedure were consistent with the results of jack- 
knifing estimates performed using 14 of the 15 fields (ie. 
an n — 1 jackknife estimate) or using 8 of the 15 fields 
(ie. an approximately ~ n/2 jackknife estimate). The 
value for the fit is 3.7 on 4 degrees of free dom. 

As noted in Sl. lAdelberger fc Steidell ()2005D performed 
a cross-correlation measurement with a similar sample of 
color-selected galaxies to compare black hole and galaxy 
masses over a large range of AGN luminosities (-20 > 
Mab(1350A) > -30) at a similar range of redshifts to our 
galaxy sample (1.5 ^ z < 3.6). That study separated the 
AGN sample into two bins of black-hole mass, obtain- 
ing galaxy-AGN cross-correlation lengths tq ~ 5.27'![_\ Iq 
h-^ cMpc for AGNs with 10^-^ < Mbb/Mq < 10^ 
and ro = 5.20l};^^ cMpc for AGNs with 10^ < 
Mbr/Mq < 10^°'^. These measurements are fairly con- 
sistent with our own measurement of ro for the galaxy- 
HLQSO cross-correlation, given the si ze of the uncertain- 
ties, and lAdelberger fc Steideil (|2005f) assume a correla- 
tion function slope of 1.6, rather than the slope of 1.5 
used in this study. As noted above, the best-fit value of 
ro varies inversely with the chosen slope for the range 
of separations our measurements include, and this effect 
likely accounts for the slight discrepancy between these 
two estimates. 

4.2. Comparison to Galaxy-Galaxy Clustering 

The strength of the clustering signal corresponds to 
the mass scale of the HLQSO-host halos, which we are 
interested in comparing to the average halo mass scale of 
non-active galaxies. As such, the relative strengths of the 
galaxy-galaxy (GG) and galaxy-HLQSO (GQ) clustering 
reveal the relative mass scales of their respective host ha- 
los, and therefore illuminate any halo-mass requirements 
for the formation of HLQSOs at z ~ 2.7. 
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Fig. 5.— The ratio E°ti='/EP™'^ for the galaxy-QSO cluster- 
ing (GQ; top) and galaxy-galaxy clustering (GG; bottom). In 

each panel, the red diamonds denote Tj°^^ /Y^y'^^'^ , while the black 
crosses denote S°'=7SP'''='^. The error bars are from Poisson uncer- 
tainties. The dashed black line is the fit to Sg'-'^/Eg'^^'^ and defines 
the angular selection function; the fun ctional form is a power-law 
f or th e GQ selection function ( i|4. 1.111 and linear for the GG case 




Projected Distance (h ' cMpc) 



Fig. 6. — Estimate of the reduced galaxy-HLQSO (red) and 
galaxy-galaxy (black) correlation functions Wp{R) for those galax- 
ies closer than 1500 km s~^ from the HLQSO (or fiducial galaxy) 
redshift. Wp is the excess probability of a galaxy appearing at a 
projected comoving separation R from the HLQSO line of sight, 
as compared to predicted galax y nu mber counts determined by 
the re dshift selection function ( qS.ll l and angular selection func- 
tion ( i|4.1.1l l. Solid cur ves a re fits to the best-matched Multi- 
Dark halo populations ( ^4.31 1. which imply a galaxy halo mass 
log(Mh_gai/MQ) = 11.9±0.1 (see ^Ot and a HLQSO halo mass 

log(Mh!QSo/M0) = 12.3±0.5. 

Our estimate of the galaxy-galaxy correlation function 
is based on the same technique as our galaxy-HLQSO 
estimates but is modified by centering on each galaxy, 
rather than the HLQSO, in turn. In addition, we restrict 
our GG analysis to those galaxies at redshifts Zgai > 2.25 
so that the GG autocorrelation function probes a similar 
redshift range to that of the GQ cross-correlation; the 
909 galaxies with z > 2.25 have a median redshift z™!*^ = 

2.63, while the median HLQSO redshift is = 2.66. 

For each galaxy in our sample, we consider the number 
density (per unit solid angle) of galaxies as a function of 
projected distance from our fiducial galaxy, separating 



between redshift-associated galaxies (those within 1500 
km s^^ and in the same field as the fiducial galaxy) and 
non-associated galaxies (those outside the velocity range 
or in a different field) . We then integrate over the redshift 
selection function fi )3.ip to find the expected number of 
galaxies in each interval, from which we can define an 
angular correlation function for each interval (by analogy 
to Eq. [n]): 



VGG 

'(;,prcd 



VGG 
^O,prod 



[1 



„GG 



(14) 



where again the v subscript denotes quantities cor- 
responding to the redshift-associated sample, and the 
subscript denotes those corresponding to the non- 
associated sample. The non-associated sample is not ex- 
pected to cluster about the arbitrary line of sight defined 
by the position of our fiducial galaxy, so we interpret the 
quantity I]obs,o(-R)/5^prcd,o(-R) as an estimate of the rela- 
tive completeness of our angular sampling. We found the 
completeness (ie. the angular selection function of the 
galaxy-galaxy pairs) of our sample to be well-described 
by a linear model in R with negative slope: So ^ aR + b. 
This shape reflects the fact that we are able to mea- 
sure the power on small scales for essentially all galaxies, 
while we can see the maximum separation ~ 2i?niax only 
for the small fraction of galaxies at the edge of our fields 
(and even then we see only the subset of pairs that lie 
entirely within the field). The best- fit parameter values 
for the GG angular selection function are a = —0.160, 
b = 1.04 with R in cMpc (Fig. [5]). 

As in the case of the GQ cross-correlation function, we 
then fit a model to that is a combination of the 

underlying clustering signal described by w^^(i?) and 
the selection function described by The combined 

model is given by Eq. 1151 



^GG 



(i?) 



Eprcd 
V 



(aR+b) 



„GG \ T 



R 



i?2 



(15) 

Fitting this model to the measured values of S^^(i?fc), 
we find the best-fit galaxy autocorrelation length to be 
r^'^ — (6.0±0.5)ft-~^ Mpc, again fixing the slope 7 — 1.5. 
In this case, the errors in S^'^ cannot be considered Pois- 
sonian because each galaxy is counted in several pairs and 
the galaxy counts are correlated between bins. However, 
the 15 HLQSO fields each provide an independent es- 
timate of S^^, and the error used for the x^-fitting is 
based on the scaled scatter among these values. The 
quoted error on r^'~^ is the la uncertainty from the same 
bootstrap and jackknife procedures described for r^"^ 
fi )4.1.1[) . The value for the fit is 9.3 on 5 degrees 
of freedom. 

This autocorrelatio n length is sign i ficantly larger 
than that found by lAdelberger et alJ (|2005bl ) [tq = 
(4.0±0.6)/i~i cMpc at z = 2.9], despite both studies 
relying on a similarly-selected set of galaxies. How- 
ever, this study is restricted to the spectroscopically- 
observed galaxies, which have a higher mean luminos- 
ity than the galaxies in the photometric sample used 
in lAdelberger et al.i ()2005bi) . and are more compara- 
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ble to the higher-luminosity sub-sample of galaxies used 
in that paper, for which the authors estimated tq = 
(5.2±0.6)/i~^ cMpc. In addition, the much larger set 
of spectroscopic redshifts used in our sample allows us 
to characterize the redshift selection function with much 
greater accuracy, as well as to restrict our analysis to 
those galaxies associated with the HLQSOs in three- 
dimensional space. For example, an injudicious choice of 
zq in Eq. [10] would lower the estimated cross-correlation 
length, either by failing to count HLQSO-associated 
galaxies (zq < 1500 km s~^) or by diluting the cluster- 
ing signal by the inclusion of the voids adjacent to the 
HLQSO in redshift space {zq > 1500 km s~^). 

4.3. Estimate of Halo Mass 

The measured clustering of the galaxies in our sample 
is primarily useful in its connection to the mass scale of 
the galaxy host halos. Because the clustering strength 
of dark matter halos is a function of halo mass, we can 
invert this relation to obtain the halo mass for a popu- 
lation of objects with a given autocorrelation length. In 
practice, we perform this inversion numerically, finding 
the population of simulated halos (for which the mass is 
known) that match the clustering strength of our galaxy 
sample. 

Using halo catalogs froin the MultiDark MDRl sim- 
ulation (IPrada et al.l 120111: a ccessed via the MultiDark 
Database of Ricbc et al.ll20lTI ). we measured the correla- 
tion length ro as a function of min imum ha lo mass My, us- 
ing the Landy-Szalay estimator (jLandv fc Szalavill993| ) 
and assuming the same power- law slope (7 = 1.5) used 
in our fit to the galaxy autocorrelation function. The 
correlation function was measured for halo populations 
of differing masses by varying the minimum Mh in steps 
of 0.05 dex; the correlation lengths for a subset of the 
halo samples are listed in Table [3l A power-law correla- 
tion function is a poor fit to the halo clustering at small 
scales due to the effect of halo exclusion; therefore, we 
restricted our fit to pairs with separations 1 < d/{h~^ 
Mpc) < 5, a range that avoids the halo-exclusion zone 
while still closely matching the range of projected dis- 
tances in our observed sample. In this manner, we find 
that our galaxy sample is most consistent with having 
a minimum halo mass log(Mh/M0)>11.7±O.l, the fit to 
which is displayed in Fig. [6] The halos in this mass range 
have a median mass log(Mh/M0) = 11.9±0.1. The sta- 
tistical error in the mass estimate is entirely due to the 
propagated error in the autocorrelation function, as the 
uncertainty in the autocorrelation function among simu- 
lated halos is negligible by comparison. 

In addition to matching clustering strengths, we can 
also attempt to match the abundances of observed galax- 
ies and simulated halos. Although our spectroscopic 
sample of galaxies is incomplete, we can compare this 
ha lo population t o the galaxy luminosity function (GLF) 
of iReddv et al.l (|2008[ ) , which corrects for incomplete- 
ness in both the spectroscopic and photometric sam- 
ples. Luminosity functions are measured separately for 
galaxies with 1.9 < z < 2.7 and 2.7 < z < 3.4, while 
our sample straddles these two redshift intervals, but 
the GLF evolves very little over this redshift range, and 
the pre dictions of either model are quite similar. Us- 
ing the 'Schechter" (1976") GLF parameters listed in Ta- 
ble 7 of .Reddv et al.. (,2008. ) . and taking our magnitude 



TABLE 3 
Clustering Properties of 
Simulated Halos 



A/T 1 m m 1 1 m 

^ > 


ACF ro'* 
(/i-i Mpc) 


XCF ro'' 
(fe-i Mpc) 


11.50 


5.1 


5.6 


11.60 


5.6 


5.9 


11.70 


6.1 


6.1 


11.80 


6.7 


6.3 


11.90 


7.1 


6.8 


12.00 


7.8 


6.8 


12.10 


8.6 


7.2 


12.20 


9.6 


7.6 


12.30 


10.5 


7.8 


12.40 


11.6 


7.9 


12.50 


12.8 


8.4 


12.60 


14.6 


8.6 


12.70 


16.3 


9.0 


12.80 


18.6 


9.1 



Halo autocorrelation length (compare 
to ^ (6.0±0.5)/i-i cMpc) 
^ Cross-correlation length with halos 
of mass log(Mjj/MQ)>11.7 (compare to 

r^'^ = (7.3±1.3)/i-i cMpc) 



limit n < 25.5 to correspond to Mab(1700A) < -19.9 
at z ~ 2.7, the Reddy et al. models predict a galaxy 
number density 0gai = 2.4 — 7.0 x lO^^ft-^Mpc"^ (in- 
cluding the la limits on (f)*). The munber density of 
log(Mh/MQ)>11.7 halos in the MuhiDark MDRl simu- 
lation is 4>sim = 4.4 X lO^^h^Mpc^^ , entirely consistent 
with the measured value of 4>gai- 

Taking the population of log(Mh/MQ)>11.7 halos to 
represent the host halos of the galaxies in our sample, we 
then estimate the mass of the HLQSO hosts by finding 
the population of simulated halos whose cross-correlation 
with the representative galaxy halos is equal to our mea- 
sured galaxy-HLQSO cross-correlation function. Again 
varying the minimum Mh (of fiducial HLQSO hosts) in 
0.05 dex increments, and using a Landy-Szalay variant 
for a cross-correlation with 7 = 1.5, we find that a 
HLQSO host halo mass log(Mh/MQ)>12.1±0.5 (a me- 
dian halo mass of log(Mh/M0) = 12.3±0.5) is most con- 
sistent with our galaxy-HLQSO cross-correlation mea- 
surement. The error in the mass is due to the error 
on both Tq and the propagated error on the galaxy 
host halo mass, since the strength of the cross-correlation 
function depends on the mass of both the HLQSO-host 
and galaxy-host halo populations. The fit to the corre- 
sponding simulated cross-correlation function is shown in 
Fig. El The MultiDark halos of log(Mh/MQ)>12.1±0.5 
have an autocorrelation length of 6— 15 cMpc, which 
we consider to be an estimate of the HLQSO auto- 
correlation length, and such halos have an abundance 
= (0.17 - 5.9) X IQ-^/i^Mpc-^ at z-2.5 in the sim- 
ulation. 

4.4. Dependence on Simulation Cosmology 

The MultiDark suite of simulations used cos- 
mological parameters based on the WMAP 5-year 
results, {nm,^A,cr3M = {0.27,0.73,0.82,0.70}, 
which are consistent with the most recen t 
WMAP 7-year results from iLarson et all ()2011l) : 
{0.276±0.029,0.734±0.029,0.801±0.030,0.710±0.025}. 
In comparison, the older (and widely-utilized) Millen- 
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nium simulation (jSpringel et al.|[2005l) used cosmological 
parameters based on the WMAP 1-year results, 
{rj„,rJA,(j8,/i} = {0.25,0.75,0.9,0.73}. We here consider 
how such a variation in these cosmological parameters 
affects the halo-matching process employed in this study. 

The parameters flm and fjg both affect halo abun- 
dances, and thus affect the halo bias and the mapping 
from clustering strengths to halo masses. iZehavi et al.l 
(|2011l ) conduct an HOD analysis on a large sample of 
galaxies and find that varying the matter density over the 
range 0.25 < ^Im < 0.3 produces only a ^2% variation in 
their clustering measurements, which is quite small com- 
pared to the statistical uncertainty in our measurements. 

However, the amplitude of the linear dark-matter fluc- 
tuations, (Tg, is tied to the clustering in a more pro- 
nounced and complicated manner. The clustering of 
galaxies in linear theory is given by the galaxy bias and 
the dark-matter clustering: ^gg(-^^) = b'^{M)^uM- De- 
creasing as decreases the value of ^dm but also greatly 
decreases the number density of high-mass halos, which 
causes the bias at a given halo mass b{M) to increase. 
For high-mass halos, the overall effect is to increase the 
clustering strength at a given mass when erg is decreased, 
which suggests that mapping halo masses to clustering 
strengths using the Millennium simulation would result 
in a shift toward larger halo masses. Repeating our halo- 
matching analysis on Millennium halo catalogs, we find 
that the best-matched halo population of galaxies has 
a minimum halo mass log(Mh,Min/M0)>12.O±O.l and a 
median mass log(Mh Miii/M0)=12.2±O.l in that simula- 
tion. The discrepancy between the two simulations is 
~3x the statistical uncertainty in the galaxy halo mass 
measurements, confirming that the clustering of these 
massive halos is quite sensitive to the chosen cosmologi- 
cal parameters. 

4.5. Relative Abundances of Galaxy-Host and 
HLQSO-Host Halos 

The mass scales of the host halos for galaxies and 
HLQSOs map to halo abundances, as described in H4.3i 
A galaxy host halo abundance of 0sim,ii.7 — 4.4 x 
10~^ft.^Mpc~^ (the MultiDark Simulation abundance of 
halos with log(Mh/M0)>11.7) and a HLQSO host halo 
abundance of <^sim,i2.i = 1-2 x 10~^h^Mpc~^ suggests 
that halos massive enough to host a HLQSO are only 
^4x less abundant than those massive enough to host 
the average galaxy in our sample; the fact that far fewer 
than one quarter of the galaxies in our sample host a 
HLQSO is a strong constraint on the duty cycle of these 
objects. However, the precise value of the HLQSO duty 
cycle depends on the number density of HLQSOs, which 
in turn depends on the choice of QSO population. 

All of our HLQSOs have luminosities at rest-frame 
1450A of log(z/L^/L0)-14, or an absolute magnitude 
M(1450A) -300 This is brighter than the luminosity 
range for which large-sample statistics ar e available in 
surveys such as SDSS and SLAQ [see e.g. iCroom et al.l 
120091 whose Mg(2;=2) is comparable to M(1450A)], but 

^ This criterion may not be satisfied for the gravitationally lensed 
object Q0142-10, and it is possible that other QSOs in our sample 
are also lensed. However, we regard it as unlikely that significant 
lensing has remained undetected in these well-studied objects, so 
the rest-frame luminosities quoted here are assumed to be accurate. 



we can obtain an estimate of the z ~ 2.7 quasar lumi- 
nosity function (QLF) b y extrapolatin g the r esults of the 
highest redshift bins of ICroom et al.l ()2009D to slightly 
higher redshifts and luminosities; in this way we roughly 
estimate the number density of M(1450A) > —30 QSOs 
to be (j)Qso ^ I0~^'^h^ Mpc""^. Integrating this den- 
sity over the total comoving volume between redshifts 
2 < z < 3 predicts ~25 QSOs in this luminosity range 
over the entire sky, suggesting that a large fraction of the 
comparably bright QSOs at these redshifts are already in 
our sample. 

Given this number density, we can extract the duty 
cycle of HLQSOs from the ratio 0QSo/0sim,i2.i — 
10^^ — 10^^, defining the duty cycle as the fraction 
of halos massive enough to host a hyperluminous QSO 
[log(L/L0)>14] that actually do host such a QSO. This 
extreme rarity with respect to the number of potential 
host halos indicates that the formation of the HLQSO 
must rely on a correspondingly rare event occurring on 
scales much smaller than those probed by our analy- 
sis, perhaps related to an extremely atypical merger or 
galaxy interaction scenario. 

4.6. Black Hole Mass vs. Halo Mass 

It is interesting to compare the host halo masses of 
the HLQSOs to the minimum black hole (BH) masses 
allowed by their luminosities under the assumption of 
Eddington-limited accretion; we will refer to this mini- 
mum mass as Mbh- The minimum BH masses for each 
HLQSO are listed in Table [2] We calculate the value of 
Mbh directly from L1450 (the value of vLi, at a rest-frame 
wavelength of 1450A): 

o 1 in~5 f -^1450 \ , r \ 

Mbh = -r-^ = 3.1 x 10 — — A/q . (16) 

47rGmpC \ Lq J 

We use Li45o in place of the bolometric luminosity 
Lboi in order to avoid the additional uncertainty in the 
bolometric correcti on. The true bolometric correc tion is 
likely to be small: Ne mmen k, BrothertonI (j2010D use a 
thin accretion disc model to predict a correction factor 
Lboi/Li450 ^ 3, which i s consistent with the empirica l 
correction estimated by iNetzer fc TrakhtenbrotI ()2007f ) 
for Lsinn and scale d by the L5100/L1450 relationship of 
INetzer et al.l (|2007f) . However, there is substantial scat- 
ter in these corrections, and it may be expected that 
Lboi/Li45o approaches unity for QSOs selected by the 
most extreme rest-UV luminosities, so L1450 is a useful 
estimate (and likely a lower limit) on Lboi- 

The HLQSOs in our sample span a range of '^Sx in 
Mbh (with the possible exception of Q0142-10), and it 
is notable that the field (HS1549+1919) with the largest 
value of Mbh is also associated with the largest redshift 
overdensity in the galaxy distribution (see column A^i5oo 
in Tabled]). However, there is no clear relation between 
iVisoo and Mbh among the other fields, and our galaxy 
samples are not large enough to comment on the varia- 
tion of Mbh with halo mass. 

The median value of Mbh for our sample is 
log(MBH/MQ) ~ 9.7. This indicates that the HLQSO 
host DM halos are only ^300-2000x more massive than 
their associated supermassive BHs, even assuming accre- 
tion at the Eddington limit. The relationship between 
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BH mass and DM halo mass is uncertain even at z ~ 
(compare e.g. [F errarese 200j, iBooth fc Schavel[2010l and 
iKormendv fc Benderll2011l ). but these BHs he weU above 
the predictions of the Mbh-c or Mbh-'^c relations for any 
reasonable mapping of the DM halo mass to the bulge 
velocity dispersion a or circular velocity Uc, as demon- 
strated below. 

Thr ee such mappings are considered in iFerrares^ 
()2002| ) , in which the halo virial velocity related to 

the circular velocity by considering — Wvir (a zeroth- 
order approximation), Vc = l.Swvir (based on o bserva- 
tiona l constraints on DM halo mass profiles by iSeliakI 
I200 2D. or Vc/v^i-c given by a function of halo mass ex- 
tractcd from the N-body simulations of iBullock et al.l 
((200 li) . These three different assumptions predict BH 
masses of log(MBH/M0) = 7.0, 8.4, and 7.5, respectively 
for a halo of mass log(Mh/M0) = 12.3, corresponding to 
Mdm/Mbh =i 2 X 10^ 8 X 10^, and 6 x 10"*. In any of 
these cases, the minimum BH masses for the HLQSOs 
in our sample are 1-2 orders of magnitude higher than 
the predictions of the low-redshift associations, implying 
that the host halos must "catch up" with the BHs in 
order to fall on the established relations by z ~ 0. 

Though estimates of BH masses at high redshift are 
highly uncertain, as are the stellar masses of their host 
galaxies, this result agrees qualitatively with several ob- 
servational studies that find BH host galaxies at high 
redshift (l<z<4)ofa given stellar mass have sys- 
tema tically higher BH masses than in the local universe 
(e.g. iPeng et al.ll2006l: iDecarh et al.l(20ToHMerloni et al.l 
[2010: IGreene et al.ll2010D : IBooth fc Schavd (I20ia (20Tll) 
describe an interpretation of this evolution in terms of 
the compactness of DM halos, which are more tightly 
bound at high redshifts. These studies generally find a 
smaller deviation from the z ~ relations than is present 
in our sample, but the extreme luminosities of the HLQ- 
SOs in our sample force us to probe the highest-mass 
end of the BH mass distribution, so our measurements 
are very sensitive to the scatter in the Mh— Mbh relation 
as well as evolution in the mean. 

If we consider the possibility that the dynamical mass 
discussed in [J5] includes matter external to the HLQSO 
host halo at z ~ 2.7, which may merge into a single 
more massive halo by z ~ 0, we can calculate where such 
a halo would fall in the Mh— Mbh relations of Ferrarese 
(2002). Taking a halo mass log(Mh/MQ) = 13, the above 
prescriptions predict BH masses of log(MBH/M0) — 8.3, 
9.6 and 8.7, respectively, which lie much closer to the 
range of BH masses seen in our sample. However, it 
seems clear that the extremely high BH masses indicate 
that the HLQSOs are atypical (with respect to the gen- 
eral population of QSOs) at the smallest scales. 

5. GROUP-SIZED HLQSO ENVIRONMENTS 

In addition to the properties of the HLQSO host halos 
themselves, it is interesting to consider the type of larger 
environment these hyperluminous objects inhabit. The 
spatial scale of the galaxy overdensities occupied by the 
HLQSOs in our sample is -0.5/1"! ^^^p^ (^^ig. H -200 
pkpc), and the peculiar velocity scale of the composite 
overdensity is (Tv,pec — 200 km s~! after accounting for 
our measurement errors (Fig. [T]). The relatively compact 
nature of the overdensity suggests that it may represent 
a virialized structure (discussed below), in which case 



the inferred size and velocity scales can be combined to 
provide a crude estimate of the mass scale of the over- 
density. The virial mass estimator can be expressed in 
terms of the 3D velocity dispersion (v^), a characteristic 
radius R, the gravitational constant G, and a constant 
a ~ 1 that depends on the geometry of the system: 

Mdyn = ai?(v2)/G. (17) 

If we approximate our group as a sphere of uniform 
density, we have a — 5/3. We can also take (v^) = 
{vD + {vl) + {vD = 30-2 ~ 3 X (200 km s-^f and R ~ 200 
pkpc from the scales above, from which we find that the 
average HLQSO overdensity is associated with a total 
mass log(Mdyn/MQ) ~ 13 - the approximate mass scale 
of a small galaxy group, and consistent with the HLQSO 
host halo mass derived from the clustering analysis in 

Because of the crude nature of this estimate, we con- 
sidered several checks to determine whether these over- 
densities are indeed consistent with virialized groups. 
If the galaxies around the HLQSOs are in virial equi- 
librium, their spatial extent should roughly match the 
virial radius r2oo of a log(Mh/M0) ~ 13 halo, where 

^200 = (3Mgip/8007rpcrit)^^'^- In fact, the virial radius 
for this mass scale is approximately r2oo — 200 pkpc 
~ 0.5/i~! cMpc - this close match to the observed over- 
density scale suggests that the HLQSO-associated galaxy 
overdensities are indeed virialized. 

Finally, we can estimate the number of galaxies asso- 
ciated with each HLQSO. From our smoothed selection 
function ( §3.1|) . we find that the approximate number 
density of spectroscopically-observed galaxies at z — 2.7 
is (/"spec = 1-3 X 10~^/i^Mpc~^, while the number den- 
sity of log(Mh/Mo)>11.7 halos in the MuhiDark sim- 
ulation is 0sim = 4.4 X 10~^h'^Mpc^^ (which is also 
the abundance of galaxies predicted by the GLF of 
iReddv et al.l I2008D . Under the assumption that each 
of these log(Mh/M0)>11.7 dark matter halos host a 
galaxy of comparable luminosity to those in our sam- 
ple (i i4.3|) . this implies that our spectroscopic sample is 
~30% complete. We find a total of 15 galaxies in our 
sample that are within 1500 km s~! and 0.5h~^ pro- 
jected cMpc of a HLQSO; taking our completeness into 
account, we expect that there are another ~ 35 galax- 
ies remaning unobserved in this volume. On average, 
therefore, each of the 15 HLQSOs in our sample has —3 
other log(Mh/M0) ~ 12 galaxies within 200 pkpc, again 
suggesting a group-sized environment. 

Note that we use the term "environment" here to con- 
note a region that may or may not correspond to the 
host halo of the HLQSO. The mass we derive here is 
slightly larger than the average HLQSO host halo mass 
of log(Mh,Qso/Mo) = 12.3±0.5 derived from our clus- 
tering analysis ( ij4.3p . and the galaxies associated with 
the HLQSO overdensity extend to greater projected radii 
than the —130 pkpc virial radius of such a halo (Fig. 
[4|). The discrepancy in the mass estimate may be due 
to larger-than-assumed errors in the HLQSO redshifts, 
as noted in §3.31 overestimation of the galaxy velocity 
dispersion would infiate the dynamical mass estimate of 
the system. However, it may also be that the HLQSO 
host and its galaxy neighbors are subhalos within a 
larger structure corresponding to our measured dynami- 
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cal mass. 

In addition, wc note that the velocity scale of 200—300 
km s~^ and the overdensity of galaxies in such an envi- 
ronment are extremely conducive to mergers and dissi- 
pative interactions among galaxies. We suggest that the 
results of this study are thus strong evidence that the 
fueling of these HLQSOs is associated with merger ac- 
tivity, with the caveat that our sample of HLQSOs are 
extreme outliers in the QSO luminosity distribution, and 
thus may be formed and sustained by rather different 
mechanisms than the average QSO at these redshifts. 

6. SUMMARY 

We have used a large sample of galaxy redshifts to in- 
vestigate the environments of 15 hyperluminous QSOs 
(HLQSOs) in the redshift range l.h < z < 2.9. Our 
galaxy sample includes 1558 spectroscopic redshifts be- 
tween 2; = 1.5 — 3.6 from the KBSS-we use the galaxies 
far from the HLQSOs to characterize our redshift selec- 
tion function in much greater detail than is possible with 
purely photometric samples. Furthermore, all the red- 
shifts in our sample are projected within ~3' of one of 
the HLQSOs, which allows us to describe the HLQSO 
environments on sub-Mpc scales. The principal conclu- 
sions of this work are given here: 

1. The HLQSOs are associated with a (5 ~ 7 over- 
density in redshift when considered on scales of 
^hhr^ Mpc. The overdensity has a velocity scale 
of CTv, pec — 200 km s^^ after subracting the effect 
of redshift errors, and a projected scale of i? ~ 200 
pkpc. When stacked at the redshifts of the HLQ- 
SOs, the combined galaxy distribution shows no 
peaks of similar significance, and stacking on ran- 
dom galaxy redshifts shows that the HLQSOs are 
correlated with much more significant small-scale 
overdensities than the average galaxy in our sam- 
ple. 

2. Careful treatment of the HLQSO redshifts is es- 
sential in order to accurately determine which 
galaxies are associated with the HLQSOs in three- 
dimensional space. When available, we used a 
combination of low-ionization broad lines, narrow 
emission lines, and the onset of the Ly-a forest in 
the HLQSO spectra themselves in conjunction with 
narrow Ly-a at small angular separations from the 
HLQSOs to obtain HLQSO redshifts offset by hun- 
dreds or thousands of km s^^ from their previously 
published values. The velocity scale of the observed 
overdensity, which is smaller than the measured off- 
set for any one of these HLQSOs, demonstrates 
both the accuracy of our redshifts and the inad- 
equacy of common techniques for estimating the 
redshifts of these hyperluminous objects. 

3. The best-fit autocorrelation function for the subset 
of galaxies in our sample with z > 2.25 (zmed — 
2.63) has a correlation length rf/^' = (6.0±0.5)/i-i 
cMpc. Comparison to dark-matter halo catalogs 
from the MultiDark simulation suggests that the 
galaxies in our sample have a minimum halo mass 
of log(Mh/M0)>11.7±O.l and a median halo mass 
of log(Mh,med/M0) = 11.9±0.1. 



4. The best-fit galaxy-HLQSO correlation func- 
tion for our sample has a correlation length 
rf^ = (7.3 ± 1.3)/i-^ cMpc. By measur- 
ing the clustering between dark matter ha- 

los of varioiis masses and those halos having 
masses log(Mh/MQ)>11.7, we find that the cross- 
correlation between log(Mh/M0)>11.7 halos and 
log(Mh/M0)>12.1 halos most closely matches our 

observed value of r^*^. We therefore deduce 
that each HLQSO in our sample inhabits a dark 
matter halo with mass log(Mh/M0)>12.1±O.5, 
which corresponds to a median halo mass of 
log(Mh,med/M0) = 12.3±0.5. The number density 
of these halos exceeds the number density of HLQ- 
SOs by a factor ^10*^ - 10^. 

5. The HLQSO luminosities imply minimum masses 
log(MBH/M0) > 9.7, suggesting a BH-DM mass ra- 
tio Mdm/Mbh ^ 300—2000 for a dark matter mass 
log(MDM/nisun) ~ 12.3—13. Such a small ratio in- 
dicates that the HLQSOs are significantly overmas- 
sive with respect to the Mbh^ Mh relation at z ~ 0, 
and appear overmassive with respect to equivalent 
relations at higher redshift (thoiigh black hole mass 
estimates are quite uncertain at high redshifts). 

6. The HLQSOs in our sample are associated 
with group-sized environments with total mass 
log(Mgrp/M0)'-^13. This conclusion follows from 
a dynamical mass estimate from the peculiar ve- 
locities and projected scale of the galaxy overden- 
sity, and is consistent with the virial radius and 
galaxy counts expected for such a group. The pe- 
culiar velocities and overdensities associated with 
groups strongly indicates that these HLQSOs in- 
habit environments where mergers and dissipative 
interactions are common. 

In conclusion, the results of this paper demonstrate 
that the host halos of HLQSOs are not rare, so the 
scarcity of these objects is likely due to an extremely 
improbable small-scale phenomenon that produces HLQ- 
SOs. Such a phenomenon could be related to an atypical 
galaxy interaction geometry or similar scenario: the over- 
dense environment with small relative velocities would 
increase the probability of such an event, but an unusual 
merger configuration is likely required to generate such 
large black hole masses and QSO luminosities. 
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